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We consider the Poisson-Nernst-Planck system which is well-accepted for describing 

Q-t dilute electrolytes as well as transport of charged species in homogeneous environ- 

<D ■ 

ments. Here, we study these equations in porous media whose electric permittivities 
show a contrast compared to the electric permittivity of the electrolyte phase. Our 
main result is the derivation of convenient low-dimensional equations, that is, of ef- 
fective macroscopic porous media Poisson-Nernst-Planck equations, which reliably 
describe ionic transport. The contrast in the electric permittivities between liquid 
and solid phase and the heterogeneity of the porous medium induce strongly oscillat- 
ing electric potentials (fields). In order to account for this special physical scenario, 
we introduce a modified asymptotic multiple-scale expansion which takes advantage 
of the nonlinearly coupled structure of the ionic transport equations. This allows for 
a systematic upscaling resulting in a new effective porous medium formulation which 
shows a new transport term on the macroscale. Solvability of all arising equations 
is rigorously verified. This emergence of a new transport term indicates promis- 
ing physical insights into the influence of the microscale material properties on the 
macroscale. Hence, systematic upscaling strategies provide a source and a prospec- 
tive tool to capitalize intrinsic scale effects for scientific, engineering, and industrial 
applications. 
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I. INTRODUCTION 

The Poisson-Nernst-Planck (PNP) equations can be applied in many different physical 
contexts such as modeling of ionic transport, e.g. batteries, supercapacitors, fuel cells, and 
capacitive desalination devices. Especially the fields of electrokinetics and electrohydrody- 
namics gained an increasing interest in recent years. Actual research aims to take advantage 
of scale effects in micro- and nano-fluidic devices for industrial applications and for the cre- 
ation of chip-like devices ("lab on a chip"). Such devices can perform separation, mixing, 
and chemical analysis tasks. It is also possible to design electrokinetic pumps.- 

The study of geometric effects on the scale of cell membranes, muscles, and neurons by 
means of PNP equations finds currently a lot of attention in biology and medcine.- 1 ? The 
essential goal is to better understand how calcium ions, i.e. Ca 2+ -ions, move in voltage- 
dependent calcium channels for example. These channels are a group of voltage-gated ion 
channels which can be found in muscles, glial cells, and neurons. Recent research attempts 
to mimic such biological ion channels with synthetically built channels.- For example by 
modifying channel geometry and surface charge one tries to better understand the effect of 
rectification. Rectificiation can be descriptively explained by the comparison of ionic flux 
with an electric current through a pn-diode. One usually studies rectfication factors (ratio 
of forward current to reverse current) in this context, see la for example. 

This broad range of applications in heterogeneous environments strongly rely on models 
which reliably and systematically account for effects of the microscale on the macroscale. A 
very common approach for deriving effective macroscopic equations is volume averaging.-~- 
Unfortunately, it is still unclear how to systematically treat nonlinear terms by this intuitive 
method. A technically slightly more involved approach is the homogenization method^"— 
which provides a reliable and systematic alternative under the assumption of a periodic pore 
distribution. 

The general importance and the strong demand of properly upscaled equations in en- 
gineering, and design as well as optimization of scientific and industrial devices call for 
mathematical tools that rely on well-established principles for multiscale problems. Here, 
we want to systematically extend the widely accepted PNP system from the free space case 
towards solid-electrolyte composites showing a high contrast between its electric permittiv- 
ities. For this purpose, we consider the full PNP equations with the help of a modified 
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asymptotic two-scale expansion. This new approach accounts for the nonlinear and coupled 
structure of the system, see Theorem IIII.21 As our main result, we derive the following new 
effective macroscopic porous medium PNP equations, that is, 

pd t u r - pAu r + div {W{t, x) Vug) - div (z r u5MVug) = in Q T , 
-div (e Vug) =p(uj- ug) in Q T , 

which is valid under local (pore level) thermodynamic equilibrium and for arbitrary Debye 
lengths Xo > 0. The parameter p denotes the porosity and W, M, and e° are effective trans- 
port tensors defined by the upscaling subsequently performed. The variable Uq represents 
effective macroscopic quantities such as the concentration of positively charged ions n + for 
r = 1, the density of negatively charged ions n~ for r = 2, and the electrostatic potential $ 
for r = 3. 

All equations appearing during the upscaling are rigorously justified by well-posedness 
criteria. In particular, Lemma |TTL4]( Section [TTT|) guarantees the solvability of the new system 
(TjQ) which shows a new term D r that accounts for a dominant influence of the oscillating 
electric potential on the concentrations. We emphasize that this new term emerges as a 
result of an adapted asymptotic multi-scale expansion introduced in order to account for 
the heterogeneity induced by the porous medium and by the electric permittivities. That 
means, in the classical asymptotic expansion 

u s (x) = uq(x, x/s) + sii\(x, x/s) + s 2 m 2 (x, x/s) + . . . , (2) 

we assume a special separation between the micro- and macro-scale in the terms u\ and 
«2 in difference to related literature^— and classical homogenization theory,— 1 ^ see ( 12^1) 
in Theorem IIII.2I below. We point out that the nonlinear character of the Nernst-Planck 
equations leads with (F2]) to ill-posed reference cell problems without additional assumptions. 
Hence, we suppose that the reference cells, which define the micro-geometry of the porous 
medium, are in thermodynamic equilibrium. This guarantees then well-posedness. This 
solvability issue might also be the reason why the upscaling of the PNP equations was 
mainly restricted to thin double layer type approximations^- 1 ^ or linearized formulations in 
the context of Onsager reciprocal relations^ 1 ^ so far.— ^ Physically, such situations occur 
when the electric permittivities between the electrolyte and the solid material are far apart, 
see Section ITAl 
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The article is organized in the following way: The dominant oscillating behavior of the 
electric potential is motivated by the contrast in the electric permittivities between the solid 
and the liquid phase in Section II Al where also a related effective media theory is discussed. 
A historical overview of closely related upscaling results is given in Section II BL In Section 
HTl we state elementary results and introduce necessary notation. The main results follow in 
Section IHH Finally, we prove all results in Section IIVI 

A. Physical motivation: Dielectric permittivities of solids and liquids 

In this section, we state the physical setting which leads to strongly oscillating electric 
potentials in composites such as a porous medium permeated by a dilute electrolyte. A 
related example where such oscillations are well-known is the electric field over a material 
with strongly heterogeneous conductivities. This is also one of the classical fields of effective 
media theory^ and homogenization theory*^ where one often assumes a periodic represen- 
tation of the heterogeneities for simplicity. In fact, the high-frequency electric permittivity 
and the low-frequency electric conductivity are formally equivalent because of the equiva- 
lence in the governing equations. However, the situation for the PNP system here is slightly 
different since we have to deal with a non-linearly coupled system of equations. As explained 
previously, we account for this difference by a non-standard asymptotic expansion that fac- 
tors the strong influence of the electric potential in. Moreover, the equivalence between 
permittivity and conductivity implies that their mathmatical computation is equivalent. 

As in the case of conductivity, the strength of the oscillations can be controlled by the 
distance between the different electric permittivities for our problem here. Since we study 
a dilute electrolyte, we can expect an electric permittivity of the liquid phase to be around 
80 at room temperature and a frequency under 1kHz (of course, this also depends on the 
electrolyte employed). For the solid phase, we can expect an electric permittivity between 
2 and 5, i.e., paper 3, alumina 4.5, teflon 2.1, porcelain 5.1, and plexiglas 2.8. But in many 
fields, a systematic derivation of effective media quantities such as the electric permittivity 
is still lacking. 20 Our subsequently derived equations reliably define such an effective electric 
permittivity for periodic porous media. We emphasize the importance of characterizing 
porous materials with respect to dielectric properties in microelectronics.— Moreover, we 
believe that a systematic and reliable upscaling of such complex composites using geometric 
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and material properties together with experimental validation gives promising perspectives 
for new scientific, technological, and industrial applications. 

B. Review of related upscaling results of the PNP equations 

We briefly give a shortened historical overview, see Table HJ and point out differences of 
contributions mainly based on the homogenization method. Table [T]examplifies how different 
upscaled formulations arise upon different assumptions initially (i.e., on the macroscale) 
made. 

TABLE I. Effective PNP equations: Various different upscaled formulations arise under different 
assumptions. We abbreviate by "PB" the " Poisson-Boltzmann equations", by "PNP" the "PNP 
equations" , and by " Onsager" the " Onsager reciprocal relations" . 





Article 


PB 


linearized PNP 


PNP 


Onsager 


Jackson ct al. 
1963^ 


• given periodic electric 
potential or PB cqua- 









• stationary case 

• effective diffu- 
sion equation, no 
homogenization 

Gross ct al. • no homogenization 

1968^ 



M. H. Holmes • generalization of Jackson 

1990^ ct al. 1963 including surface 

charge 

• homogenization in Slot- 
boom variables 



• surface charge 



Moync 


et 


al. 






• similar to Ray et al. 2011: 








and 






Kind of thin double layer 




Moync 


et 


al. 






limit (PB in reference cells) 




2004^4 










• surface charge 




Looker 


et 


al. 


• 


local equilibrium 






2006^ 






• 
• 
• 


local clcctroncutrality 
fluid flow 
surface charge 




s 


Allaire 


et 


al. 


• 


local equilibrium 




• rigorous 


2010^, 


based 


• 


bulk clcctroncutrality 




with two-scale 


on Looker et 


al. 


• 


fluid flow 




convergence 


2006^£ 






• 


surface charge 






Schmuck 2011 


15 






• local equilibrium 




Schmuck 2012 


2r> 






• without surface charge 




Ray 


et 


al. 






• Dcbyc length as homoge- 




2012^ 










nization parameter 

• without surface charge 

• fluid flow 
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Based on the stationary Nernst-Planck equations and a given periodic electric potential, 



22 



In 



23 



a model for 



the authors derive an effective diffusion equation for ion densities in 
charged membranes separating dilute aqueous salt solutions is studied. The authors give 
expressions for Onsager's reciprocal relations^ 1 ^ without assuming small Debye-lengths nor 



22 



a Debye-Hiickel linearizaion of the Poisson-Boltzmann equation. Based on ideas from 
Holmes^ performs an asymptotic two-scale expansion for a PNP formulation rewritten in 
Slotboom variables. This work gives a very interesting approach for nonlinear diffusion 
in charged polymers. The results seem to be closely related to upscaling of electrostatic 
potentials </>o(jz) om y depending on the micro-scale y (which is related to the use of periodic 
potentials in I22I ) 

In ns 



a macroscopic electrokinetic formulation describing electro-osmotic flow and elec- 
trophoretic motion in periodic porous media is obtained by the classical multiple-scale expan- 
sion method. They also perform a kind of thin double layer approximation in the reference 
cells by Poisson-Boltzmann equations. This approximation is well-known and frequently ap- 
plied in electro-chemistry— The same authors^ apply the periodic homogenization theory 
to upscale the Nernst-Planck equations in the medical and biological context of cartilage. 
13j, the fundamental Onsager reciprocal relations^ 1 ^ together with positive definite- 



In 



ness of corresponding upscaled tensors are derived under the assumption of local thermo- 
dynamic equilibrium. Their starting point is a linearized^ iV-component electrolyte in a 
dilute Newtonian fluid (small Reynolds number) flowing through a periodic porous medium. 
The authors also consider the physically interesting case of surface charge. However, the 
influence of the surface charge is not obvious at the end. The impact of such charges on 
the macroscopic level is of major scientific and engineering interest in microfluidics^ - — and 
neurobiology^. 

Very recently, Allaire et al.— put the physically relevant derivations from [l3| into the 



rigorous framework of the two-scale convergence. The main purpose is again the verifica- 



tion of Onsager's reciprocal relations as in [13|. Derivations of such relations require the 
assumption of local thermodynamic equilibrium, a linearized PNP system,— and an elec- 
troneutrality assumption in the bulk which is physically closely related to a thin double layer 
approximation. 

Inll4j. the authors perform a singular limit with respect to the dimensionless Debye length. 



A weighted Debye length, i.e., A a , and the use of A as the homogenization parameter has 
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the meaning of upscaling the PNP system parallel to taking special (a > arbitrary) thin 
double layer limits of the system. Espeically for a — 1, this is an interesting problem since 
the thin double layer approximation is a widely used simplification as already mentioned 
above. 



II. NOTATION AND PRELIMINARIES 



The following classical exposition recalls central definitions and results from l9Hl ll. For 
the microscopic variable y := -, we obtain the following relation for gradients applied to 
functions ip 8 (x) := %j) (x, |), i.e., 

V^ 8 (x) = V^O, y) + -V y tP(x, y) . (3) 

Homogeneous Neumann problems for Poisson equations for example require the use of 
the quotient space Wjj(y) := H}(Y)fR. This space of equivalence classes is defined with 
respect to the relation, 

u ~ v <^ u — t> is a constant, Vit, v E Ha(Y) . (4) 

For notational brevity, we do not introduce additional notation for an element of this equiva- 
lence class. A representative element of this equivalence class can be chosen by the following 
mean zero condition, that means, 

W S (Y) := | it g H}{Y) I M Y (u) = 0} , (5) 

where, 

My(u):=^rJ^u(y)dy. (6) 

Lemma II. 1 The following quantitiy, 

\W\\ Wf = \\Vu\\ L2{Y) \/u,v eW^Y) , (7) 
defines a norm on W$(Y). Moreover, the dual space (W$(Y))' can be identified by the set, 

{Fe(H}(Y))\F(c) = VcGM} , (8) 
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with, 

Definition II.2 Let c,C 6l be such that < c < C and let D C E^. We call N x N 
matrices A = Wij}^ j <N G (L°°(U)) NxN strongly elliptic, if for any u G and a.e. in -D 



(A(x)w, m) > c |w| 2 



(10) 



In our analysis we mainly have to deal with A = {iy} Ki j <JV which obvously satisfies the 
conditions of Definition III. 21 

Theorem II. 3 Let A be a strongly elliptic matrix with Y -periodic coefficients and f G 
(Wjj(y))'. Then the problem, 

Find u G WtfK) such that 

(11) 

(AVu, Vv) Y = (f, v) Wv G W${Y) . 
has a unique solution. Moreover, 

\\u\\ Wi0r) < -||/|| wy)) • (12) 

Remark II. 4 Since Theorem \II. 3\ makes a uniqueness statement, we consider in this case 
the space W$(Y) defined in |3p. We apply this convention in the whole article. 

We frequently use the following space, 

V(flr) := £°°(]0, T[; L 2 (Q)) n L 2 (]0, T[; H\Q)) . (13) 

A. Review of the classical PNP equations 



Before we come to the main results in this article, we briefly recall basics about the PNP 
system. In view of computational convenience (block matrix solvers, e.g. 34), notational 
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clearity and compactness, and a non-linear (i.e., non-symmetric) extension of the classical 
Onsager relations,—^ which classically only hold in the linearized case,— ^ motivate us to 
recall the PNP equations from Il5l written by field vectors u := [n + ,n~, $] as, 



A§u = I(u) . 

u(0, x) = h . 



u 



u 



for (t,x) e n T ■= [n T ,n T ,n T y 

infi := 

on := [T l T ,T l T ,T 



[u)V n uY := Sij {u)Viu j = 



T 

on T r rr 



on 1 j, :- 



i y 



(14) 



[r^r^r' 1 ' 



Tj 5 



D 



for % = 1, 2, 3 



where x := [x, x, x]' with x G f2 C K corresponds to the coordinate field for each component 
of the field vector u and t := [t, t, t]' with t e]0, T[ for any T e R+ is the accordingly defined 



time field. The variables u 1 



n 



u 



and u 3 = $ represent the concentration of 



positive ions, the density of negative ions, and the induced electric potential, respectively. 
We further use the convention := ]0, T[xQ. The notation accounts for the fact that the 
components of the field vector u are defined in different domains of the porous medium later 
on, i.e., either in the whole domain Q or only in the electrolyte phase Q s . We further denote 
A§u := div (S(u)Vu) with S(u) := (u)} 1<ij<d for the field indices 1 < i,j < 3, 

1 < k, I < N 

the coordinate indices 1 < k, I < N and Sj fc ^(u) = Sij(u)5ki with 5m the Kronecker symbol, 
V„, := n • V with n the normal vector pointing outward of Q, '■= [T 
:= T l T U and T^, T L T for t = r,l are correspondingly defined, and 



£,r£,r£]' with 



D t := 



d t 
d t 




{s ii (u)} 1 < i)i < 3 := 



div := Idiv 



div 
div 
div 



1 n+ 
1 -n 
A 2 



V :=IV := 



V 
V 
V 



u:= [n + ,n-^]' , I(u) := [0,0, 
gi := [nf,nY, (pi]' , g r := [n+ , n~ , <p r ] ' , h := [h^h^O]' . 



n + — n 



(15) 



r 1 represents the Dirichlet (for t = D) and Neumann (for t = AT) boundary surrounding the 
porous medium Q, i.e. dQ = T D U P^. Hence, the first equation dH])i is equivalent to the 
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following classical system, 



d t n + = div (Vn + + n + V$) , 
d t n~ = div (Vn _ — n~V$) , 
-A 2 A$ = n + - n~ , 



(16) 



which can be interpreted as a gradient flow of the following free energy, 



F = U -TS 




^u^logu'-^ + ^^u^-A^V^I 2 dx. (17) 



i i 



We recall that ( fTTj) builds the basis of dilute solution theory which accounts for thermo- 
dynamic quantities such as entropy S formed by the first summand in the integral (fTTj) . 
The remaining integrands such as energy density of interactions (second term) and energy 
density of the electric field (third term) constitute to the internal energy U. We note that 
from the energy (fTTj) we can obtain the chemical potential of the ion densities u 1 = n + and 
u 2 = rT by taking the first variation with respect to u 1 and u 2 , respectively. 

An interesting question is whether the minimization of the free energy ( fTTj) by a gradient 
flow also follows the physically relevant path far from thermodynamic equilibrium. In which 
physical sense does the flow with respect to the Wasserstein distance*^ - — provide optimality? 

III. THE MICROSCOPIC POROUS MEDIUM FORMULATION AND 
MAIN RESULTS 

The study in this article relies on the system ( fTij) reformulated for periodic porous media. 
A scaling parameter s is defined as the ratio between the microscopic length scale I and the 
macroscopic size L of the porous medium, i.e. s := j- <C 1. It is assumed that s scales the 
periodicity of the reference cell Y C M. N which defines the micro-geometry. In this reference 
cell, we denote the fluid (liquid) region by Y s C Y such that its complement is the solid 
phase. After periodically covering the domain Q by such cells Y, we denote the resulting 
macroscopic domain of the periodic union of the subsets Y s by Q s and its complement by 
B s := Q\Q S . Hence, the perforated domain B s represents the solid phase and fl s the liquid 
phase, see Figure [TJ 
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Reference cell Y 



Periodic covering by scaled cells Y 





@ 




r^i r^i r rj i r^i 














@ 




r^l r^l r rj i r^l 







Homogenous approximation 




FIG. 1. Left: Strongly heterogeneous material as a periodic covering of reference cells Y := [0,^]^. 
Top, middle: The reference cell Y represents a characteristic mean pore geometry. Right: The 
"homogenization limit" s := — Y stands here for the leading order approximation of non- 
standard two-scale expansions. 



Under these considerations, the material tensor § from (fbil) ? depends now on s too, i.e., 



<i,j<3 



1 n+ 
1 -n- 

e(x/s) 



>kl 



for 1< k, I < N . 



(18) 



where e(x) := X 2 Xn s (x) + cxxb s (x) with the classical dimensionless Debye length A := ^f- 
of the PNP system ( fT6|) . a = is the dimensionless dielectric permittivity where e m and 
€f are the dielectric permittivities of the solid and liquid phase, respectively. In ( fTSl one 
recognizes that also the physical quantities like concentrations n+, n~ and electric potential 
depend on the scaling parameter s. Hence, the problem f fT4|) reads now in the periodic 
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setting as follows, 

B t u s - div (§ s (u s )Vu s ) = I(u s ) , in Cl s T := [Q S T , Q S T , Q T }' , 

u s (0,x) = h s , in n s : = [ft s ,ft s ,fi]', 

§ s (u s )V n u s = 0, onT?, 

u s = g L , onTy, t = l,r, (19) 

(§ s (u s )V n u s )f l = 0, on I*, fori = 1,2, 

(§ s (u s )V n u s ) f 3 |^= (§ s (u,)V n u s ) f 3 | fif , on/', 

^|n.= , on , 

where fj := [^1,^2,^3]' f° r i = 1,2,3 and 5^ is the Kronecker delta and Ij* := dQ? \ 
U j the solid-electrolyte interface. From ( !T9|) it follows that the flux with respect to 
u is in general not differentiable. This motivates to study f lT9|) in the sense of weak solutions. 
Moreover, the main difficulty and difference of this work is the nonlinear structure which 
prevents the material tensor § s (u s ) to be a strongly elliptic operator. For convenience, we 
rewrite (jl9~j) for [n+, n~, $ s ] = [u*, u^, ujj] in the classical form here, that is, 

d t n+ = div (Vn+ + n+ V$ s ) in tt s T , 

d t n~ = div (VnJ - n s "V$ s ) in Q S T , (20) 
— div (e(a;/s)V < l ) s ) = nj" — in Q T , 

where the corresponding ionic fluxes are defined by, 

j+ : V//.; • //JV'I'.. 
j; := Vn; - n;V$ s , 

and the boundary and initial conditions are, 



(21) 



h] n 7 = hl inn*, 

n t =9i n~ = 9l on r T > ( 22 ) 
j+-n = j;-n = oiiT?, 

where 1 = l,r. Our main result relies on the assumption of local thermodynamic equilibrium, 
which is widely used and generally accepted.- 
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Definition III.l (Local thermodynamic equilibrium) We say that the reference cells Y are 
in local thermodynamic equilibrium if and only if for each x/s = y of the same reference cell 
Y it holds that 



where /ig denotes a constant value of the chemical potential of positive (r = 1 ) and negative 
(r = 2) ion densities. Hence, the locally constant potential /Zq only assumes different values 
in different reference cells. 

In order to simplify the subsequent derivations, we make the following 

Assumption (ABC): (Academic Boundary Conditions) We assume for the ion densities 
u r f orr = 1 ; 2 homogeneous Dirichlet boundary conditions on dQ and no-flux boundary con- 
ditions on the solid- electrolyte interface I s := dQ s \ |r D U T N } inside the porous medium Q. 
For the electric potential u^ = $ s; we assume homogeneous Neumann boundary conditions 
on the boundary of the porous medium, i.e., d£l, and continuity, i.e., &\ n3 = &\ BS ; at the 
solid- electrolyte interface I s := dQ s \ {T D U T N } as well as continuity of the corresponding 
fluxes, i.e., e(x/s)V„$ s | Qs = e(x/s)V n $ s | Bs on I s . 

These considerations allow us to state our main result which consists of passing to the 
limit s — > 0. This limit has the physical meaning of homogeneously mixing the solid and 
the liquid phase in the porous material under constant volume fractions, see Figure [TJ Such 
a homogeneous description of a porous medium is a good approximation if the medium is 
very heterogeneous, i.e., s = jr is very small. 

Theorem III. 2 We assume that the reference cells Y are in local thermodynamic equi- 
librium. The boundary conditions satisfy the Assumption (ABC). Then, the solution 
u s (t,x) := [uj(t, x),ul(t, x),u^(t, x)] of problem admits the following formal asymptotic 
expansions, 



/io := logu r (t,x) + z r u 3 (t, x) = const. , 



(23) 




8 ]T C k (t, X, X/S) ^ + S 2 C M (t, X, X/S)UI + ... ft 

k=l k k,l=l 



"or r = 1,2, 



(24) 
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where C k (;;y) e V(n T ,W t (Y s )) solves e k (y) e W tt (Y) sote gfl, C r "(v,y) e 

l^(fi T , VFjj(y s )) so/ves ( fyg| ), ( 3kl (y) e Wjt(Y) so/ves ( [7^ ), and u is a solution of the following 
upscaled system, 



pd t u r - pAu r + div (D r (t, sc) Vug) - div (^u^MVug) = zn fi T 



-div (e Vu3) =p(uj- ug) 



in f2 T , 



(25) 



where p := |Y S | / |y| is the porosity. The tensor W(t,x) := {D kl (t,x)} 1 



<k,KN 



related to 



diffusion, the tensor M := {M^} 1<fcKAr related to electro-convection, and the effective per- 
meability tensor €° := {^} 1<S . KJV are defined by, 

N 



1 

1^1 ./ys 



M 



1 

IF, ./v 



i=i 



6; 



5^{^-%9 % .e 3fc (y)} dy, 

i=i 

- / £e(y) (^-(Sy^Cy)) 

1 >/y 7=1 



(26) 



The second order correctors ( rkl only exist for positive densities Uq. 

We note that for the classical asymptotic two-scale expansions*^ an upscaling is per- 
where also error erstimates between the microscopic periodic formulation and 



formed in 



25 



the upscaled equations are derived. 



15 



Remark III. 3 (1) Theorem UII.^ is an extension of the two-scale convergence results in 
by the non-classical asymptotic expansion j2^ x ■ We note that the expansions [Efy are only 
formal because convergence of such series is a priori not guaranteed and possible boundary 
layers are neglected. 

(2) The effect of the upscaling in the above theorem can be best seen in the change of the 
material tensor (Q2j), which reads for the new system ( T131) as follows 



{s' ihjl {x,Uo)} 1 



<i,j<3 



forl<k,l<N. (27) 



pS M -D^ + z r ulM M 
p5 ki -D£, - z r u 2 M M 
e° kl (x/s) 

A comparison of pT7| ) with l[T8\) clearly motivates the use of the term "material tensor" in 
the context of porous or composite media. 



14 



New porous medium Poisson-Nernst-Planck equations 



(3) The effective material tensor (2l\ ) can also be considered as a generalized effective, con- 
centration dependent conductivity tensor (as in heat/ diffusion equations). 

We note that the effective macroscopic material tensor (|27|) reveals that one can cont- 
role and enhance material transport by adjusting the contrast in the electric permittivities 
between the electrolyte and the porous medium. The different upscaling ideas as well as 
different physical and mathematical assumptions play an important role for an improved un- 
derstanding of how different microscopic material properties such as pore geometries, electric 
permittivities, charge numbers of the ions, and surface charges of the porous medium influ- 
ence the transport characteristics on the macroscale. 

Our last result guarantees the existence and uniqueness of solutions for the porous media 
PNP equations fl25|) . In order to keep the presentation clear and to prevent a technical 
derivation, we restrict our considerations to purely academic boundary conditions in the 
next 

Lemma III. 4 Let Q C M. N be a domain with smooth boundary dQ and T > small enough. 
The boundary conditions satisfy Assumption (ABC). Then, the coupled system of equations, 

pd t u r - pAu r + div CD r (t, x)Vu 3 ) - div (^uJJMVu 3 ) = in Or, 

(28) 

—div (e°Vu 3 ) = p (u 1 — u 2 ) in Q T , 

together with the following initial and boundary conditions, 

u r (0,x) G L 2 {VL) f) Hq(Q) for r = 1,2, 

u r = ondttx]0,T[ for r = 1,2, (29) 
V n u 3 = ondQx]0,T[, 
has unique weak solutions u r G V(Q T ) := L°°(]0, T[; L 2 {Q)) n L 2 (]0, T[; and u 3 G 

L°°{]0,T[]H 2 {Q)). That means, u r {t,x), r = 1,2, andu 3 {t,x) solve 

p| (u r , p r ) + V (Vu r , Vy^ r ) - (D r (t, x) Vu 3 , V(^ r ) 

+ (z r u r MVu 3 , <p r ) = V</ G H^(Q) , (30) 

^ (e°Vu 3 , V<^ 3 ) = p (u 1 - u 2 , v? 3 ) V^ 3 G H\n) , 
where h\q) := {u G H\Q) \ J Q udx = 0}. 

From a rigorous point of view, Lemma rill.4l finally guarantees that the second order terms 
in the asymptotic expansion fl24l) are locally well-defined. 

15 
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IV. PROOF OF THEOREM UTOI 

We first rewrite the second term in equation (fl~9]) i with help of the property ([3]) in the 
following way 

div (§m.x)VV> s ) = * 2 div„ (S (if) s , y) Vyip s ) 

+ s' 1 {div, (S(V>„, y .lV.V.s) + div v (S(^„ y)V,V,) } 

(31) 

+ div., (§(^,y)V x ) 

:= [s- 2 5o + .s '5, + <S 2 ] ^ (x, *) , 

where V>s( x ) is an arbitrary function as in Q. For the subsequent considerations we apply 
the following notation, 

[S 2 ] N ij> s := [div(S(^,y)V^)f := [div (S(^, y)V v 1>.)] N 

■= ( div * M*> y)^i) ^ii^M^) (32) 

= div x (^V^ s 3 )fi - div^V^s^ • 
Hence, the operator [•] represents a restriction to nonlinear components given by the indices 
(i = 1, j = 3) and (i = 2, j = 3). 

We make now the formal Ansatz of the asymptotic expansion 

u s (t,x) « u (t,x, x/s) + sui(t,x, x/s) + s 2 u 2 (t,x, x/s) + . . . , (33) 
with Uj(t, x, y) for i — 0, 1, 2, . . . such that 

Ui(t, x, y) is defined for (t, x) G^andyGY := [Y s , Y s , Y)', 

(34) 

Uj(-, •, y) is Y-periodic. 



The above Ansatz is formal because there is no guarantee that the series (133]) is finite. 

(1) Problem for terms of order (D(s~ 2 ): After inserting (1331 into (fT9lb . using fl3T|) 
and (1521 gives a sequence of problems by equating terms with equal power in s, that means, 

{— S(un, y)u = in Y , 

1 } (35) 
u (-,y) Y-periodic. 

If we use definition (13T|) . then we can rewrite (!35|) i as the following equation, 

-div (S(uo, y)V,u ) = in Y , (36) 
16 
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wich is equivalent to the system, 



0{s 



-A^uJ - div^ (uj V y v? ) 



-A y u 2 + dw y (u 2 Q V y u 3 Q ) = 
-div y (k(2/)V v u§) = 



in y s . 
in Y s . 
in y . 



(37) 



One recognizes immediately that solvability must be first established for equation ( 1371) ?. 
This is immediately achieved by Theorem 111.31 Moreover, this theorem implies that Uq(x, y) 
is invariant (constant) in y G Y as a solution of (1371)?. i.e., 



ul(t,x,y) = ul(t,x) . 



(38) 



Using invariance ( 1551) in equations ( 1571b and ( 1571) 9 implies with Theorem III. 31 the additional 



mvanances 



Uo(*, x, y) = Uo(t, a;) and ug(t, ac, y) = Uo(t, ar) 



(39) 



Let us go over to the next problem in the sequence of equal power in s. 

(2) Problem for terms of order 0(s^ 1 ): (Reference cell problems) The second problem 
has the form, 



Ois- 1 ) : 



-<So(u ,y)ui = [<So(ui,y)] Ar u + «Si(uo,uo,y)uo in Y . 



Ui(-, u) y-periodic 



(40) 



Let us write (140]) in a more intuitive form by its single components, i.e., 



-A y u\ - divy (uoV y u?) = div y (uJV y ug) + div^VyuJ 



Ois- 1 ) 



+div x (uoVj/ug) + div^V^Uo + div y (uJV^u 



-A y ul - divy (ugVyU?) = -divy (ufVyuf,) + diVa-V^ 



-div^ (uqVj/U^) + div^V^Uo - div y (u§V x u; 
-divy («(2/)V J/ Ui) = div^ (n(y)V y ul) 



+div y (K(y)V x u, 



in Y & 



in y s 



in y 



(41) 



The system ( 14T|) is a linear, elliptic second order partial differential equation. Hence solv- 
ability of (|4ip follows immediately with Lax-Milgram's Theorem by starting with problem 



17 



New porous medium Poisson-Nernst-Planck equations 



[T]) 3 . The fact that u is independent of y together with the linearity of (j4"Tj) and that Sq 
only contains derivatives in y, motivates to make the following Ansatz for U!(t,x, y), that 
means, 

N 

ul(t,x,y) = -J2C'(t,x,y)d Xj ul(t,x) for r = 1,2, 

( 42 ) 

ul(t,x, y ) = -J2^(y) d ^ u U^ x ) ■ 

We use now ( )42l) and the independence of uo of y in order to rewrite ( 14T1) as a problem for 
£ Tj for r = 1, 2 and 1 < j < N as follows, 

(*, x, y)d X] nl - Y^=i d yi {<\ i (<l,.u' l) + z r u r d Xj u 3 )) 

= -dw y (z r n r V y e^y)) d Xj n 3 in ffxP, (43) 

[-div, («(y)V^(y)) + ZlA(<y)tii) = in Q x Y. 

We point out that under local thermodynamic equilibrium, that means, in each reference 
cell we have due to the induced separation of scales by the limit e — > 0, 

see Definition IIII.ll Hence, the term with the summation on the left-hand side in (T4"3"l) i 
disappears. System (1431) defines the reference cell problems for the porous media corrector 
functions £ rj for r = 1, 2 and 1 < j < N . Such correctors finally define the effective tensors 

(USD. 

Remark IV. 1 We point out that the Ansatz [JW ^ s an extension from linear homogeniza- 
tion theory and is canonically chosen to account for the problem's coupled and nonlinear 
structure. The interpretation of the Ansatz ffifyi is that oscillations in the microscopic vari- 
able of the electrostatic potential dominate the oscillations of the concentration variables. 

Lemma IV. 2 Let Uq G V(Qt)- There exists a unique solution £, 3j (y) G W^Y) for each 
1 < j < N of problem f(3\)2, i-G., 

Find & e WAY) such that 

* (44) 

a? (e>,w) = F?H := Eti {K(y)5 lv d m w)y ^ G W t (Y) , 

18 
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where af(^,w) := (n(y)V y ^, V y w) Y . With f l (y) G Wjj(Y) for 1 < I < N also the 
existence and uniqueness of a solution t; r i(t,x,y) G V(Qt, W$(Y s )) of problem [(3\)i follows, 
that means, 

For each 1 < j < N find ^(t,x,y) G V(Q T , W t (Y s )) such that 

(45) 

(^, w r ) = F[ 3 (w r ) Ww r G L 2 (Q; W^Y 8 )) , 
where for all w r = 4> r (y)ip r (x) G Hq(Qt, W$(Y)) we define, 

aI(r j ,w r ):=((V w r'(*,x,!/),V^ r )y.,V r ) n , , , 

(46) 

F?K) := - ((:,!,;, V.C'-i//). V y f) YS ,V) n ■ 
Proof. 1 The lemma is a consequence of Lax-Milgram's thoerem. 

Step 1: Problem (14411 : The assumptions of Lax-Milgrams are easily verified since n(y) is 
a strongly elliptic matrix, see Definition \II.B . 

Step 2: Problem ([45]): a) Continuity: Forv r , w r = <p r {y)4j r {x) G H^(Q,W t {Y s )) we can 
estimate the bilinear form &[ by, 



|al(v r ,w r )| = |((V w v r ,V w p ) y .,^ r ) n | < IMI L 2 (fWs(ys)) H r \\ Wf{ Ys) WWmn) > 
and hence continuity follows. 

b) Coercivity: Forv r , w r = (fi r (y)ip r (x) G Hi (Q, W$(Y s )) we derive a lower bound by, 



(47) 



% (v , v j 



((V^V,v% s ,l) Q = ||V,v 



r||2 

lL 2 (n ; L 2 (y s )) 



(48) 



c) Ff* is linear and continuous: For w r = (fi r (y)?p r (x) G i/p(f2, 14-^(1^)) we estimate F^ by, 
|F?K)| < |((^ V^,V^)y-^) 

<c||v^ 3 ' 

and continuity follows with £ 3j G W^(y) obtained in Step 1. 



(49) 



ft" 1 ||i2(y s ) || L 2(y s ) ||U || L 2(^) 



Il 2 (h) 



(3) Problem for terms of order 0(1): We collect all the terms which do not contain 
any factor e. As a result we end up with the following equation 



0(1): 



-«S (u , y) u 2 = «Si(u , ui, y)m 

+5 2 (uo, ui, u 2 , y)u + I(u ) + D t u 
U2(-,y) ^-periodic, 



(50) 
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where Si := Si(u ,y) + [<S (ui, y)]^ and S 2 := <S 2 (u ,y) + [^(u^y)]^ + [<S (ui,y)] . In 



order to study solvability of ( 1501) it is an advantage to write (15011 explicitly for each physical 
quantity. For r = 1, 2 we have 



Oil) : ^ 



— = dlVy (ZrVL^VyUl) 

+ {diy r (zrX^Vyuf) + div y (z r v^V x ul) 
+div x V y \f 1 + div y V x VL[ + div^ (z r UiV y uf)j 
+ {A x u r + div x (z r u r V x ul) + div x (^^V^u 3 ,) 
+div y (z r u[V x ul) + divy (zrU^VyUo)} — d t u r Q in Q s x 
-div y (/c(y)Vj,u|) = div x («(y)V 2/ uf) + div^ (/c(y)Va;uf) 

+div 2; (/«(y)V a; ug) + (uj - ug)xy= in V . 



(51) 



Since a solvability constraint implies the effective equation for the macroscopic quantities, 
we first achieve the well-posedness for the system (jSTl) . 



Lemma IV. 3 The problem [51\) o. i.e., 

Find u|(- , • , y) G W$(Y) such that 

(52) 

4 (ul w 3 ) = F 3 (w 3 ) Vu> 3 G W B (Y) , 

has a unique solution u 3 G H^-^Qt, W$(Y)) where we define, 

a 3 (4,w 3 ) := (K(y)\7y4,\7 y w 3 ) Y , 

F 3 (w 3 ) := (div x . (K(y)VyU 3 ) ,w 3 ) y - (K(y)V x u 3 ,V y w 3 ) Y (53) 
+ (div^ (K(y)V x u 3 ) ,w 3 ) y + (u£ - ug, w 3 ) ys . 

For r = 1, 2 equation /[5l\) i . that means the following problem, 
Find uJ^-, G Wj|(F s ) s,uc ^ 

(54) 

a£ (u5, w r ) = FI>(w r ) Vw r G ^(F 5 ) , 
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has a unique solution u^-, -,y) G ^(F 5 ), where we for z\ = 1, z 2 = — 1 define, 



4(u> r ) := (V tf u5,V 



yW )y s 



F r 2 (w r ) := - (z r u r V y nl, VyW r ) Ys + (div x (z r u r V y ul) ,w r ), 

- (zruJVxU?, V y w r ) Y3 + (div x Vyu[,w r ) Ys - (d t u r ,w r ) Y3 (55) 

- {V X U[, VyW r )y s - (ZrU[VyUl,VyW r ) Ys + (A^Uq, W/")y s 



+ (divx (z r vT V x u r ) ,w r ) Y3 - (z r u[V x \il, V y w r ) 



ys ■ 



Proof. 2 Step 1: Theorem UI. S\ and Lemma \II.l\ immediately provide existence and unique- 
ness of a solution G IV^F), if we have, 

(^2, ^)(Hi(Y°)y,Hi(Y°) = • (56) 
Equation 156]) reads in terms of physical quantities as, 

(div x . (k(j/)V„u?) , l) y - (div a (K(y)V^) , l) y + (uj - ug) ys = . (57) 

27ms equation defines the upscaled formulation (f^Sj)? /or i/ie electric potential such that l[5b}) 
holds true. 

Step 2: Again, solvability follows by Theorem UI. Si after verification ofF r 2 G (Wfl(Y s ))'. .Due 
to Theorem \II. 3\ and Lemma \II. 1\ it must hold for r = 1, 2 that 

l)(H B 1 (y s ))',H s 1 (y s ) — , (58) 
which reads in explicit form as follows, 

(diV^VyUi, l)y s + (A^Q, l)y. - (^Uq, l)y, 

(59) 

+ (div x (z r u r V y u\) , l) y . + (div x (z^VyvQ , l) ys = . 



Since we use / fJgj) represents the by Lemma \III-4\ well-posed effective model, we herewith 
guarantee that 155]) holds. 



With representation ( l42l) 9 we can rewrite (J57l) in the following way, 

TV 

- ("M (*« - M!) > !)r = l y 1 ( u o - u2) • (60) 

i,j,k=l 
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Using (1551) allows us to write ( 1601) more intuitively by, 



ME^ = mK-uS, (6i) 

i,j,k=l 

such that e° := {e^Ji^AKJV is defined by 

4 : = - 14 ((*(v) (fo - t.Ar)) , i) y • (62) 

We can write down in the same way effective equations for equations (1541) by using (|42|) i . 
That means we rewrite (l59l) in the following way, 



(^(^^Uo), l) ys + (d Xi ({-Sijdy.^ "} d Xk u 3 ) , l) Y3 - (d t u r , 1] 
+ (d Xi (z r u r {5 ik - 5 tl d yj e k } d Xk u 3 ) , l) y . = . 

If we apply (I59|) we can even further simplify (l65j) to 



(63) 



(64) 



\Y S \ d t u r - \Y S \ Au r + d Xi H {Mw^fo v)} dyd Xk u 3 c 

- d Xi (z r u r J {5 ik - 'i/'^.C"" (//)} dyd Xk v% \ = . 

As a consequence of ( I64I) we obtain the effective diffusion-related tensor D r : — {D^} 1<i k <Ni 
i.e., 

j^jf^Mv^frz.y)} Vi, fc = 1, . . . , N , (65) 
and the electro-diffusion-related tensor M := {M^}^ k<N , i.e., 

M ^ := uTT f - Mw£ 3 *(?/)} ^ Vz, = 1, . . . , iV . (66) 

The definitions (165]) and (166]) finally provide the porous media approximation of the Nernst- 
Planck-Poisson equations, 

{pd t u r - pAu r + div (© r (t, x) Vug) - div (z r uj;MVug) = , 
(67) 
-div (e°Vu[j) = p(ul - u 2 ) , 

lyS I 

where p := W> is the porosity, r = 1, 2, zi = 1, and Z2 = — 1. 
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(4) Derivation of the second order correctors: In order to compute the second 
order corrector for u 3 , we use (JUJh in equation flEQh , i>e., 

AT 

-dw y {n(y)VyU 3 ) = - E d *i «y)^j9 yj ^ k d Xk nl) 

i,j,k=l 
N 

- E M^MM^M!) 

i,i,fe=l 

+ div^ (k(v) V^Uq) + (uj - Uq)xys in F . 
Inserting equation ( )6T|) into (!68|) leads to the following problem, 

<9 2 u 3 



(68) 



N 



-dw y (K(y)V y ul) = -^2e° kl - 



d 2 u 3 Q 



N 



dxudxi 

k,l=l i,j,k=l 



E K Wkj9 Vj i 3 



dxidxk 



N r) 2 n 3 N /} 2 n 3 



(69) 



' dxjdxi 



in F. 



i,j,k=l , Z = 1 

With equation f l6"$j) the right-hand side F 3 in ([SSI ? can be rewritten by, 



N 



(F2,w 3 )( H / tl (y))',iy t ,(Y) = E 



- /■ 



E / /<y)4A (£ 3 '-^)™ 3 ^ 



(70) 



The same arguments as those for ([4"2]) suggest to look for a function u 3 , of the following form, 



N 



ul(t,x,y):=J2( 3kl (y) ^ 



fci=i 



dx k dxi 



(71) 



where C 3fei is the solution of, 

-div, ( K ( y ) v,c 3 -) = -4 - e ^ m^m*^ 31 ) - E <v) 6 kAi (£ 3i - Vi) in y • 

(72) 



N 



»J=1 



3=1 



Lemma IV. 4 There exists a unique ( 3kl for 1 < k,l < N that solves equation (72) in the 
weak sense, that means, ( 3kl is a unique solution of the following problem, 



Find ( 3kl G W t (Y) such that 



4(C 



3ki w 3\ 



? 3 fe! i 



W" 



Ww 3 e W t {Y). 



(73) 
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where we define. 



3 (C 3fci ,^ 3 ) := (K(y)V tf C 3 ",VX) 



Y 
N 



iV 

-E(#AP"i» 3 )r' 



(74) 



j'=i 



Proof. 3 T/ie existence and uniqueness is an immediate consequence of Theorem ill. 3\) and 
Lemma (II. 



The same considerations as those for the derivation of (j72]) can be applied to (155]) (or to 
( l5Tj) in the context of its classical formulation). With ( 171]) and (J4"2"]l we can rewrite ( )5T|) as, 



fe,z=i fe,(=i fe,i=i 

j,k,l=l i,j,k=l J j,l=i 



- V 1*5,. 9 f^^i + ^a.a £3 

j,k,l=l 



dxidx k 



(75) 



i,j,k,l=l K 1 i,j,k=l 

i,j,k=l J i,j,k,l=l 



7,Z=1 k J ' 



dxjdxi 



in F s . 



We make now a corresponding Ansatz for the functions u£ with the indices r = 1,2 under 
the same considerations as those for (ITT]) , i.e., 



iV 



uS(*,z,jO= ^C fci (M,2/K- (76) 
fc,i=i 

After inserting definition (ITS"]) into ( 175]) we obtain an equation for the second order corrector 
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functions ( rkt , that means, 



N 



N 



fcj=i 



k,l=l 



' N N 

P A— 1 » 1 



J = l 



J^uJ_ 

dx k dxi 



N 



N 



iV 

+ E (<W fe Kc* - i)) 



E (^ 3i - f 



k,l=l L 3=1 
TV 

E ^ (Mjw^-^C 8 ")) 



9ug Suq 



2 r UQ 



d 2 ug 
dxhdx 



N 



k,l,j=l N J 



»A7 



dx k dxi 



(77) 



In order to guarantee solvability of equation ([77]) . we assume that Uq > for Uq(0, x) > 77 > 0. 
This can be obtained by special test function techniques as applied in^i to prove non- 
negativity of solutions. Existence and uniqueness of the corrector functions defined by 
equations ( I72l) and ( 1771) is achieved in the following 



Lemma IV. 5 Let \x r Q G V(Qt), u o(t,x) > for all (t,x) G Qt and 1 < r < N , and assume 
that the reference cells Y are in local thermodynamic equilibrium, see Definition UII.1\) . 
Then, there exists a unique ( rkl for 1 < k, I < N that solves equation in the weak sense. 
That means, ( Tkl is a unique solution of the following problem, 



Find C kl (t,x,y) G V(Qt, Wfi(Y s )) for 1 < k,l < N such that 



a£(e*W) 



Vw r G Hq($1, W$(Y s )) . 



(78) 
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where we define forw r := <p r (x)ip r (y) , 



a5(C«,«; r ) := ((uS«(y)v s C rw ,v y ^)y.^ r ) n , 



N 



N 



K) := ( ( DL - E 4A - w) - E ^ (M^f fc ) , ■ 

j=l ij=l 

AT 



+ -M H -^^.p-jft) 



i=i 



N 



+ E ^ (Mwf fc - 1)) 



r Z r 



N 



d 2 < 



i=i 



E^ (M«(^-W)/ 

»,J'=1 



<9uq <9uq 



r? 2 u 3 



dxkdx 



d 2 u 3 
dx k dxi 



(79) 



Proof. 4 PFe appiy Lax-Milgram's theorem in three steps: 
a) Continuity: For v(x,y),w(x,y) G ifg (Q; W$(Y S )) we have 



\a 2 {v,w) 



K(y)VyV,V y w) Ysn <C\\V y v\\ Hl ^. L2{Ys)) \\Wyw\\ Hl{n . L2{Ys)) , (80) 



which proves continuity since \\^yv\\ L2 (Y«) defines a norm on W$(Y) and k G M(c,C,Y). 
b) Coercivity: Forv(x,y) G Hq(Q; W$(Y s )) we can estimate 



a£(u,u)|= (u5«(?/)V w u, V w u) y . jn > c ||V w u||^ 1(n . La(y . 



(81) 



c) F^' is linear and continuous: In order to keep the clear representation we only study 
the two critical terms, i.e., 

(/) =- - D.^)v. (gg) .r) 

W^e need to s/iow too* 1(7)1 , |(i7)| < C || 

w lli : / 1 (QT;W / t |('K s )) ■ The second term (II) can be esti- 
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mated for w r (x,y) = (j) r (x)ip r (y) G Hq(Q; W$(Y s )) by, 



ICTI<C ||e 3 i Hl(n;L2(ys)) + ||V y C 3 -IL^)llV, 



<9 2 u 3 



dxkdxi 



(83) 



u 



L2(Q) 



£/ie regularity available for £, £ and Ug ; the continuity of the second term (II) follows. 
We estimate the first term (I) by, 



(84) 



I (-01 < C \\C k \\ H i^ Wi ( Y ^) ( || V ^ 3i |L3(y 3 ) 



+ 1J ll^ 7 l/ < / ,r |lL2(y S ) || U o||H2 (f7 ) ll U olli/l(n) 11^ HfTl(fi) ' 

which implies the desired continuity. 



V. PROOF OF LEMMA lHL4t WELL-POSEDNESS OF THE UPSCALED 
MODEL 

For each v r G V(Qt), t — 1, 2, consider the linear system 



p<9 t u r - pAu r = -div (D r (i, x) Vv 3 ) + div (z r v r MVv 3 ) in VL 7 



(85) 



in 



-div (e°Vv 3 ) = p (v 1 - v 2 ) 

where B' (t, x) indicates that its t and x dependence originate from v r for r = 1,2. The 
choice of v r G V(Qt) guarantees that the right-hand side in (l85lb is in L 2 (Q). Hence, there 
exists a unique solution u r G L 2 (0, T; ifg(f2)), <%u r G L 2 (0, T; L 2 (fi)) by standard parabolic 
theory. 

In the same way, take v r G V(£lr), r — 1, 2, and let u r solve, 

p<9 4 u r - pAu r = -div (© r (t, x) Vv 3 ) + div (z r v r MVv 3 ) in Q T 
-div (e°Vv 3 ) = p (v 1 - v 2 ) in tt T 



86) 



After subtracting (186]) from (185]) . we obtain the following equation for u r := u r — u r , 



])0,r ._ 



- J°' r , for r = 1,2, i.e., 

pd t u r - pAu r = -div (D°' r (t, x) Vv 3 ) - div ( D°' r (t, i)Vv 3 

+z r div (v r MVv 3 + v r MVv 3 
-div (e° Vv 3 ) = p (v 1 - v 2 ) 



in fi T , (87) 
in Qt ■ 
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The test function u r , r = 1,2, in (157)1 together with Sobolev inqualities and L p -interpolation 
estimates induce the following inequality, 



p d 

-— I llu 

2dt 



in 2 



;2||2 



a 



in 2 



Vu 



2 ll 2 



<C7(a,e,^ir)(||v 1 |r + 
which is a consequence of the following estimates, 



+ e(||Vv x | 



|Vv 2 | 



(88) 



^0,rv7— 3 



Vw\ Vu' 



v 3 l 



< c m 

< Cllv 3 



H2 



max 

l<k,l<N 



D 



,1/2 



I// 2 



|Vv r ' 



L3 
,1/2 



|Vu r 



|Vu r 



(89) 



< C(a,e) v 



2 , a 



li/ 2 



v r r + e Vv r z + -||Vu 



and 



Vv 3 , Vu r 



< C(a) ||v r " 2 



t (||Av 3 || ||Vv 3 ||) + ^||Vu 



< C(«) ||v r ||^ ||Vv 3 || (llu 1 )! 2 



^2 1 1 2 



- Vu 2 



(90) 



It is immediately clear that the terms containing M can be controlled by the same bounds, 
since M is just a constant matrix. After integrating (|88|) with respect to time we get, 



P 
2 



in 2 



\L°°(0,T;L 2 (n)) 



+ U 



2 II 2 



lL°°(0,T;L 2 (n)) 



HI 2 



lL 2 (0,T;L 2 (n)) 



+ Vu 2 



lL 2 (0,T;L 2 (f2)) 



< C(e,a,v r ,v r )T (Wv 1 ] 
+ e (ll V ^l 



L°°(0,T;L 2 (n)) 

IVv 2 ""' 



+ v z 



lL°°(0,T;L 2 (n)) 



(91) 



lL 2 (0,T;L 2 (n)) 1 II v " llL 2 (0,T;L 2 (n)) / 

By choosing T and e such that C(a, v r , v r )T < p/4 and e < a/4, we obtain a unique 
solution by Banach's fixed point theorem. The statement that u 3 G L°°(]0, T[; H 2 (Q)) is an 
immediate consequence of elliptic regularity theory. 
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